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In this article, we show that eigenenergies and eigenstates of a system consisting of four one¬ 
dimensional hard-core particles with masses 6m, 2m, m, and 3m in a hard-wall box can be found 
exactly using Bethe Ansatz. The Ansatz is based on the exceptional affine reflection group Fi 
associated with the symmetries and tiling properties of an octaeube —a Platonic solid unique to 
four dimensions, with no three-dimensional analogues. We also uncover the Liouville integrability 
structure of our problem: the four integrals of motion in involution are identified as invariant 
polynomials of the finite reflection group Fa, taken as functions of the components of momenta. 


Introduction.- The relationship between exactly solv¬ 
able quantum one-dimensional multi-particle problems 
and kaleidoscopes has been long appreciated [l[ . A kalei¬ 
doscope is a system of mirrors where none of the multi¬ 
ple images of the original object is broken at the mirror 
junctions. A viewer situated inside a kaleidoscopic cav¬ 
ity has no means of distinguishing the images of objects 
from the original. For a broad class of boundary condi¬ 
tions, the eigenmodes of a kaleidoscopic cavity are repre¬ 
sented by finite superpositions of plane waves ^hey can 
be found exactly, using the method of images 3 . When 
a many-body problem reduces to a solvable kaleidoscope, 
the resulting solution is known as the (coordinate) Bethe 
Ansatz solution [^, . The list of particle systems solv¬ 

able using Bethe Ansatz is so far exhausted by: equal 
mass hard-core bosons, on an open line, on a circle or 
in between walls [^; J-interacting, equal-mass, generally 
distinguishable particles on a line and on a circle [a-lsj; 
this includes bosons which could also be in the pres¬ 
ence of one or two walls [ij; systems of two hard-core 
particles with masses m and 3m interacting with a wall, 
for both wall-m-3m and wall-3m-m spatial orders were 
briefly commented on in Q, but discarded as inferior to 
the problems with finite strength interactions. 

It has been long conjectured that no exceptional— 
specific to a given number of spatial dimensions— 
kaleidoscopes lead to solvable particle problems with re¬ 
alistic interactions [l|, |^. The search for physical real¬ 
izations was limited to systems of particles of the same 
mass, with a possible addition of immobile walls. In this 
Letter, we show that the exceptional closed kaleidoscope 
Fa induces a novel quantum integrable system: four one¬ 
dimensional quantum hard-core particles with different 
masses, 6m, 2m, m, and 3m, in a hard-wall box. The so¬ 
lution utilizes the symmetries of an octacube, a Platonic 
solid unique to four dimensions. 

Identifying the particle problem generated by the Cox- 
eter diagram F^. — Consider six hard-core (i.e. impene¬ 
trable) particles on a line, with masses mo, mi, m2, m3. 
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FIG. 1. The Fa kaleidoscope and its particle realization. 
From top to bottom: the affine Coxeter diagram Fa; the an¬ 
gles between the mirrors; the fonr particles with masses 6m, 
2m, m, and 3m between two hard walls at xq = —L and 
xs = 0; the particle-particle and particle-wall hyperplanes of 
contact in all three coordinate systems used in the main text. 


m4, and m^. Their coordinates will be denoted as Xg, xi, 
X 2 , X 3 , X 4 , and X 5 respectively. The natural coordinate 
transformation 


Xi = 



for 1 = 0, 1, ..., 5 


reduces the system to a single six-dimensional particle of 
mass A4 (to be fixed later) whose motion is constrained 
by the following set of five inequalities: yi+i jy/rnfifi > 
yijyjrnl, for z = 0, 1, ..., 4. 

Consider now the affine Coxeter diagram Fa 0 de¬ 
picted in the top line of Fig. [TJ It encodes the geometry of 
a particular simplex-shaped 4-dimensional kaleidoscope. 
The vertices label its five mirrors: the 3-faces (repre¬ 
sented by 3-dimensional simplexes) of the 4-dimensional 
simplex. The edges—the links—of the diagram (and 
their absences) encode the angles between the corre¬ 
sponding mirrors, i.e. the angles between their normals. 
When two vertices are not linked by an edge, the respec¬ 
tive mirrors are orthogonal to each other. An edge with 
no numbers above it corresponds to an angle of 7r/3. A 
number n above an edge would produce an angle of zr/n 
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between the corresponding mirrors. The meaning of the 
mark inside the leftmost vertex will be explained later. 
In our case, the non-right angles between the mirrors are 
7 r/ 3 , 7 r/ 3 , 7 r/ 4 , and 7 r/ 3 , in the order of their appearance 
on the Coxeter diagram, left to right. 

Let us try to identify the five consecutive vertices of 
the Coxeter diagram with the five hyperplanes of con¬ 
tact between neighboring particles: for 

i = 0 , 1 , ..., 4 (see the particle diagrams in Fig. [T|). 
This choice is natural: contact hyperplanes for two non¬ 
overlapping pairs of neighboring particles are indeed or¬ 
thogonal to each other. Then the four edges of the dia¬ 
gram correspond to the four triplets of consecutive par¬ 
ticles. 

For three consecutive particles, i, i + 1, and i + 2, the 
angle between the i vs * + 1 and z + 1 vs i -I- 2 contact 
hyperplanes is given by 


(i+i) {i+i) — arctan ^ 


I rrii+i [mi + mj+i + mi+ 2 ) 

mjmj +2 


(see e.g. @). Notice that these angles do not depend 
on the overall mass scale. Therefore, any constraints on 
these angles can only lead to relationships between the 
five mass ratios, mi+i/mi for z = 0, 1, 2, 3, 4 alone, with 
no overall mass scale involved. On the other hand, for 
our six-body system to lead to the F 4 kaleidoscope, the 
angles between the hyperplanes of contact must satisfy 
four equations only. 
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thus being seemingly underdetermined. Nevertheless, 
once the non-negativity of the masses (mi > 0 for 
z = 0, 1, 2, 3, 4, 5) is invoked, a single solution (up to 
an overall scale m) survives: 


mo = -boo; mi = 6 m; mi = 2 m; 
m 3 = m; m 4 = 3m; = +00 . 


In retrospect, that both the leftmost and the rightmost 
mass should be infinite could have been predicted from 
the onset: this is probably the only way to map six par¬ 
ticles on an open line to a single 4-dimensional particle 
in a closed cavity. 

Without loss of generality, we can place the infinitely 
massive particles at 

xo = -L\ X5 = 0 , 

where L is the size of the resulting hard-wall box (see 
the hyperplanes of contact in terms of the original x- 
coordinates in Fig. [T|). Now observe that m 3 = m and L 
appear to be convenient mass and length scales, and h a 
natural scale of action. From now on, we will be using 
the system of units where 

m 3 = L = h = 1 . 


Also, it will become clear from what follows that the 
choice AI = 2 mi allows us to respect the existing con¬ 
ventions on the size of the octacube. 

Figured] also shows the corresponding formulae for the 
contact hyperplanes expressed through the transformed 
coordinates z/i, z/ 2 , 2 / 3 , and 2 / 4 . These hyperplanes form 
the boundaries of the F 4 kaleidoscope we have in mind: 
-I < \/2yI ; yi < 2 / 2 ; 2/2 < ; V^ 2/3 < 2/4 ; 

2/4 < 0 . 

The coordinate transformation z = Tz-^y ■ y, with the 
rotation-inversion matrix Tzt^y given by 

-V 6 -y/2 -I -V3\ 

-a/6 V2 1 V3 

0 -2a/2 1 a/3 

0 0 -3 a /3 / 




brings the domain of our F 4 kaleidoscope to the conven¬ 
tional form BSEHil: 


Vg = such that ag ■ z > ~ 2 ^o ■ 
and 0 : 1 , 2 , 3.4 • 2 > o| , 


( 1 ) 


where y = (yi, y2, yg, Vi), and 2: = {zi, Z2, zg, 24). 

According to this convention, the kaleidoscope dTj) 
is a 4-dimensional simplex with five 3-faces (mirrors) 
defined by five inward normals rjj = OLj/\oLj\ for 
j = 0, I, 2, 3, 4, with oo = (—1, —1, 0, 0), oi = 
(0, -bl, —1, 0), 02 = (0, 0, -bl, —1), 03 = (0, 0, 0, -bl), 
04 = (+ 5 , — 5 ,—^, — 5 ) being the so-called minimal 
(or negative of maximal) root (j = 0 ) and the simple 
roots (j = L 2, 3, 4) of the corresponding finite reflection 
group F 4 [l6| . The mirrors from I to 4 pass through the 
origin; the 0th mirror (marked by a dot in the Coxeter di¬ 
agram of Fig. d]) passes through the point ( 1 / 2 , 1 / 2 , 0 , 0 ). 
The hyperplanes that define the inequalities o are also 
identified in Fig.dl In Eq.dl the subscript the subscript e 
in /Ce has the meaning of the identity element of the finite 
group F 4 , a convention whose meaning will become clear 
from what follows. Note that the “natural coordinates” 
2 : can be expressed through the particle coordinates x as 
z = Tzi^x ■ X, with 
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Periodic tiling hy consecuitive reflections .— The key 
to Bethe Ansatz solvability of the closed kaleidoscopic 
billiards lies in their ability to periodically tile the full 
space—with neither holes nor overlaps—via consecutive 
reflections about their faces @,[l3|- Consider the kalei¬ 
doscope o. A union of the figures produced by sequen¬ 
tial reflections of © about its 3-faces from the I-st to 
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the 4-th (“simple root” mirrors), 

C = yjg^pJ^g 

= such that (gao) ■ -2 > -^(ffao) • {goto) = -1, 
for all g G F 4 I , 

( 2 ) 


with 

Vg = jz such that (goto) ■ z > —^(^ao) ■ (ffCKo) = —1 
and igai, 2 , 3 ,i) ’ 2 : > o| , 

(3) 


leads to the so-called octacube, otherwise known as the 
24-cell [ 1 ^; it is the only 4-dimensional Platonic solid 
that does not have any 3-dimensional analogues. Its 
24 vertices lie at points given by all coordinate per¬ 
mutations and sign choices applied to (± 1 , 0 , 0 , 0 ) and 
(±i, ±i, ±i, ±i). Its 24 octahedron-shaped 3-faces are 
centered at all coordinate permutations and sign choices 
of (±i, ±i, 0, 0). Sequential reflections about the “sim¬ 
ple root” mirrors of the F 4 kaleidoscope form the sym¬ 
metry group of the octacube, the reflection group ^ 4 . 
(The same label is used to mark both an open kaleido¬ 
scope, where the 0 th, “minimal root” mirror is absent 
and the corresponding simple, “non-afSne” Coxeter di¬ 
agram.) The union in ([2|) runs over all elements g of 
the group F 4 . The “physical” domain m is the identity 
member of the set 

The 4-dimensional space can be periodically tiled 
by identical octacubes; the centers of the octacu- 
bical cells lie on all points with integer coordi¬ 
nates whose coordinates sum to an even number: 
{{zi, Z 2 , Z 3 , Z 4 ) : zi+ Z 2 + Z'i + Z 4 = even}. The very 
same tiling can be obtained by a consecutive reflection 
of the o simplex about all five of its mirrors, including 
the “minimal root” one. Transformations generated by 
these reflections form the symmetry group of the octacu- 
bical tiling— the affine reflection group T 4 , homonymous 
to the corresponding closed kaleidoscope. 

In order to facilitate the possible classification of en¬ 
ergy levels by symmetry, we choose the lattice vectors 
Q' 1 , 2 , 3,4 to be proportional to the four simple roots of the 
reflection group D 4 [l^, closely associated with F 4 : ai = 


(1,-1, 0,0), a 2 = (0,1,-1,0), 03 = (0,0, 1,-1), 
and 04 = (0, 0, 1, 1) . The reciprocal lattice vec¬ 
tors are thus ki = (27r, 0, 0, 0), K 2 = {2tt, 2-k, 0, 0), 
/«3 = (tt, tt, tt, —tt), and H 4 = (tt, tt, tt, tt). 

Finding the eigenenergies and eigenstates of the 
problem .— In general, finding the eigenenergies and 
eigenstates of a H-dimensional kaleidoscopic billiard 
amounts to solving a system of D nonlinear algebraic 
equations—the so-called Bethe Ansatz equations ii 
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However, in the case of hard-core interactions, 


the problem greatly simplifies (see Corollary 3.2 in [II)- 
The eigenstates acquire the form 


= /T 7 -| 7 u exp[i(pfe)z] , (4) 

V'T’ICI -gg 


where g are the members of the finite reflection group 
G(= F 4 in our case) generated by sequential reflections 
about the “simple root” mirrors of the kaleidscope in 
question; V{g) is the parity of the number of the ele¬ 
mentary reflections needed to reach g from the identity 
transformation; Vx, = Vc/|G|(= 1/576 in our case) is 
the kaleidoscope volume; Vc = |det[(ai, 02 , 03 , ...)]| 
is the unit cell volume; |G|(= 1152 for F 4 ) is the 
order, i.e. the number of elements, of the gro up G. 
The conjecture contained in Theorem 1 in [l5| indi¬ 
cates that the wavefunction (|4|) is normalized to unity: 
/.p_ dzi(iz 2 dz 3 ... |^(z)P = 1 when integrated over the 
corresponding kaleidoscope ((P) in the case of F 4 ). 
This conjecture covers not only the case of the hard-wall 
boundary conditions but also a much more general class 
of boundary conditions, associated with finite strength 
interactions. The wavevectors k are simply drawn from 
the reciprocal lattice of the corresponding tiling [l^. In 
order to prevent both double counting and the formal ap¬ 
pearance of the eigenstates that are identically zero, the 
wavevectors k must be further restricted to those lying 
inside the open kaleidoscope associated with the kalei¬ 
doscope in question, obtained by removing the “minimal 
root” mirror: T]j ■ k > 0 for j = 1, 2, 3, ... (see the end 
of Sec. 3 of [l^b 


For the case of the hard-wall boundary conditions, the 
normalization formula used above can also be proven 
directly without resorting to the conjecture [l5|- Ob¬ 
serve that the octacube @ is a unit cell of the lat¬ 
tice spanned by the lattice vectors 04 , 2 , 3,4 and that the 
“seed” wavevector k, along with all its images gk belong 
to the respective reciprocal lattice. In that case, any two 
distinct plane waves with wave vectors gk and g'k will 
be orthogonal to each other if integrated over the whole 
cell: 


{gk\gk)c = Sg',gVc , 


with 


{ki\k 2 )A= / d‘^zexp[i{k 2 - ki)z] , 
JA 


where A is an area of space. The normalization integral 
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in question can now be evaluated as 
[ d‘'z|^(-l)^«)exp[z( 5 fe)^]p 

gGG 

9 ' 9 

r-l-g'EE(-l)^«^")+^(^)( 5 rfc|fffc)^, 

9 " 9 

= ^{99'k\gk)T>^ 

9 " 9 

9 " 9 "' 

= J2{-lf^^''H9''k\k)c = 

9 " 9 " 

= Vc = 14 )|G| ; 



FIG. 2 . (color online). Exact energy spectrum for four parti¬ 
cles of mass mi = 6m, m2 = 2 m, m3 = m, and mi = 3 m in a 
hard-wall box of length L (red solid line). Weyl’s law predic¬ 
tion for the spectrum of this system (blue dashed line).The 
energies up to 2 x lO'^/i^/maL^ are shown. 


this result justifies the choice of the normalization fac¬ 
tor in the expression dH). Above, we used (— 1 )'P (99 ) = 
(_l)T(s)(_l)T(g") and (-1)2T(9) = 1. 

We are now in the position to write down the spec¬ 
trum of our system explicitly. Starting from this point 
we abandon the m3 = L = h = 1 system of units. 
The eigenenergies, E = (/i^/ 2 A 4 )fc^, and the eigen¬ 
states 0, with k — X]j=i,2.3,4; ■'^i=i,2.3,4 G and 
%=i,2,3,4 • A; > 0, are then given by 


E., 




nin2n3n4 
2 


(2n2(ni -I- n2 + na -I- ni)-!- 


n 


GrusL"^ 

1-1-713-1-714-1- runs + 771714 -I- 773774) 


^711712113714 (^) 

~ A qa 7 4 ^ 6Xp[7(5fc„j„2n37i4) ■ 

for 772 > 1 and 771 > 774 > 773 > 1, 


( 5 ) 

( 6 ) 

•a;] 


where 


k 


"1«2713714 


/ 2771 + 2772 + 773 -b 774 \ 

2772 + 773 - 1-774 ^ ^ 

773+774 L ’ 

\ -773 +774 J 


the sum in runs over all 1152 members g of 

the reflection group +4 0 (transformations gen¬ 

erated by consecutive reflections about the four 
“simple root” mirrors of the +4 kaleidoscope); 

^nin2n3?i4 (^) = ^\ det[T)2;.^a;] | '0 nin2n3n4 ifz-i^x ■ 3^) 

are the eigenstates of the system expressed through 
the particle coordinates xi, X2, X3, X4, normalized as 

/a;i=-L/a;2=xi /a:3=X2 /a:4=a:3 = ll 

(SS) Jacobian matrix of the 2; •«— a; transforma¬ 

tion. 


At high energies, the spectrum converges to the Weyl 
law prediction for the number of states with energies be- 
low a given energy E: N (E) = 

where yV{E) is the classical phase-space volume occu¬ 
pied by points with energies below E, and D is number 
of spatial dimensions (see Fig. n. 

The ground state energy is 


E, 


ground state 


= K 


ni=3, n2 = l, n3 = l, 714=2 


l^TT^h^ 

2m3L2 


( 7 ) 


The ground state wave function can be obtained from 
Eq. (in|), using 


iround state — 3 yj2 —1,113 —1,114—2 


5 
3 

V1 y 


E(8) 


Fig. O shows a particular section of the ground state den¬ 
sity distribution within the +4 simplex, along with its 
space-tiling mirror images. 

Integrals of motion. — The Bethe Ansatz integrability 
of our system can be also reinterpreted in terms of the 
Liouville integrability. To construct the three additional 
integrals of motion in involution with the Hamiltonian 
and each other, we suggest invoking the invariant poly¬ 
nomials of the group in question [20|: finite polynomi¬ 
als of coordinates w{zi, Z2, 23, . •.) that remain invariant 
under the group action. 


w{gz) = w{z) . 


( 9 ) 


Consider now an operator / that is constructed by taking 
an invariant polynomial w as a function of the momenta 
associated with the 2: coordinates: 

i = w(-iVz) ■ 
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FIG. 3. (color online). The ground state of our four-body 
system. The plot Is performed using the 2 coordinates, rele¬ 
vant to the tiling. The image is produced by first intersecting 
the 4-dimensional density distribution by the 3-plane and the 
3-sphere indicated in the upper left and upper right corners, 
respectively (the formulas, however, are given in x coordi¬ 
nates, in which the 3-sphere becomes a 3-ellipsoid). The re¬ 
sult is a 2-sphere centered at the origin. Considering now this 
2-sphere as living in the usual 3D space, we pick, in this 3D 
space, a (2D) plane that does not intersect the 2-sphere, and 
project onto it the 2-sphere’s hemisphere closest to the plane; 
this projection is what is shown in the image. The circles 
on the 2-sphere are the sections of the 3D hyperplanes of the 
original 4D space; the great circles come from the hyperplanes 
that pass through the origin, and the other circles come from 
the hyperplanes that do not. The lighter triangle in the mid¬ 
dle corresponds to a section of the “physical” 4-dimensional 
simplex to which the particle coordinates are bounded; the re¬ 
mainder of the sphere is a section of the entire 4-dimensional 
space when this space is octacubically tiled with the mirror 
images of the original “physical” simplex. The existence of 
this tiling ensures the applicability of the Bethe Ansatz @]. 


Let us show that any eigenstate (|3]) of the kaleidoscope 
associated with the group in question is at the same time 
an eigenstate of /. Indeed, 

/>(z) = 

u;(—zVs) ^(— eyiYi[i{gk)z\ = 

9 

exp[i{gk)z] = 

9 

w(k)j2 (— exp[i{gk)z] , 
a 

Q.E.D. For the reflection group F 4 , the four lowest power 


functionally independent invariant polynomials read 21| : 




{zi - Z 2 y 

+ (zi + 

Z 2 y^ 
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II 
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II 

00 

Im= 4 : 

= 

12 



( 10 ) 


Accordingly, the four fundamental integrals of motion for 
our system read 




im 

M = 1, 2, 3, 4 




( 11 ) 


with 


, , d d , d ^ d . 

Px = —) ( 12 ) 

UXi 0X2 UX3 0x4 


being the particle momenta, corresponding to the particle 
coordinates xi, X 2 , X 3 , X 4 . Above, (.. .)^ stands for the 
transpose of a matrix. Note that the first integral of 
motion is proportional to the Hamiltonian of the system: 

= UAH. 

Summary and outlook.- In this Letter, we obtain— 
using Bethe Ansatz—an exact expression for the eigenen- 
ergies and eigenfunctions of four hard-core particles with 
mass ratios 6 : 2 : 1 : 3 in a hard-wall box. The Ansatz 
is induced by a hidden symmetry of the system related 
to the symmetries of the tiling of a 4-dimensional space 
by octacubes. The exact spectrum stands in good agree¬ 
ment with the approximate Weyl’s law prediction. 

The following observation may serve as a seed for a 
longer research program. The procedure, outlined in this 
Letter, for identifying a few-body problem relevant to a 
particular Coxeter diagram [l9| is not unique to F 4 . Any 
diagram, affine or not, that does not have bifurcations 
can potentially be used to generate a solvable few-body 
hard-core problem. TheAAr_i (N identical hard-cores on 
a circle) and Cn {N identical hard-cores in a hard-wall 
box) diagrams have been already successfully introduced 
in the first years of many-body Bethe Ansatz, in and 
[ij, respectively. The l 2 {n) diagrams were explored in 
Ref. [22j| , albeit classically, still awaiting a quantum treat¬ 
ment: for each integer n a continuous one-parametric 
family of solvable three-body mass ratios (on a line with 
no boundary conditions) was obtained. A system of two 
hard-core particles with mass ratio 3:1 in a box (whose 
one-wall version was explored in Ref. @) is expected to 
exhibit an exact solution associated with the G 2 Coxeter 
diagram. It can be conjectured that and H 4 Coxeter 
diagrams lead to solvable four- and five-body problems 
on a line, respectively. 

The true challenge is posed by the Coxeter diagrams 
(affine or otherwise) with bifurcations, most notably the 
Ee, 7 ,8 series. At the moment, it is not clear if there 
are any realistic many-body problems that can be solved 
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using these symmetries. (Some examples of “unphysical” 
many-body realizations of Coxeter diagrams—e.g. those 
where a given particle can interact with an empty point 
in between the other two—are given in Ref. Q.) 


The potential empirical context for this and further 
planned explorations is one-dimensional cold gas mix¬ 
tures in optical lattices where the dynamics is governed 
by an effective (rather than physical) mass, controlled at 
will 231. 
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